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A quantitative analysis of a microscopic model for the intrinsic Josephson effect in high- 
temperature superconductors based on interlayer tunneling is presented. The pairing correlations 
in the Cu02-planes are modelled by a 2-D Hubbard-model with attractive interaction, a model 
which accounts well for some of the observed features such as the short planar coherence length. 
The stack of Hubbard planes is arranged on a torus, which is threaded by a magnetic flux. The 
current perpendicular to the planes is calculated as a function of applied flux (i. e. the phase), and 
- after careful elimination of finite-size effects due to single-particle tunneling - found to display a 
sinusoidal field dependence in accordance with interlayer Josephson tunneling. Studies of the tem- 
perature dependence of the supercurrent reveal at best a mild elevation of the Josephson transition 
temperature compared to the planar Kosterlitz-Thouless temperature. These and other results on 
the dependence of the model parameters are compared with a standard BCS evaluation. 

PACS numbers: 74.20.-z,74.20.Mn,71.10.Fd 



I. INTRODUCTION 



The high-Tc superconductors (HTSC) reveal a number of unusual properties. One unique example, important 
both from a fundamental and from an applied point of view, is the intrinsic, "microscopic" Josephson effect: a single 
crystal like Bi^Sr^C 'aC^Os (BSCCO) consists of natural stacks of thousands of Josephson junctions, with the CuOi 
layers acting as superconducting electrodes and the BiiO-s layers as insulating ("intrinsic") barriers, closely packed 
on a microscopic length scale (d w lA). This intrinsic Josephson effect was first discovered for BSCCO- [^2) and 
later confirmed in YBaCuO- || materials. From the point of view of the microscopic theory of the HTSC, this is an 
important effect for the c-axis transport and possibly the pairing theory, ft is also an interesting and promising effect 
from an applied, technological point of view, using the microscopically packed Josephson junctions as small scale, 
tunable high-frequency oscillators 

Strong, Clarke and Anderson argued in a series of papers Jsj-Q that the unusual c-axis resistivity data obtained for 
cuprate superconductors are the result of the non- Fermi- liquid nature of the in-plane (C1LO2) ground state of these 
materials. This non-Fermi-liquid nature is taken as responsible for disrupting the interplane hopping of electrons so 
strongly that single electrons effectively do not hop between the planes, giving rise to anomalous c-axis transport 
properties || and to the anomalously large superconducting transition temperatures for these materials |(|. In this 
theory, both of these effects depend on the absence of normal interplane hopping. It is therefore natural to study the 
microscopic Josephson coupling between CuC^-planes, which exhibits this characteristic. 

A theoretical description of the intrinsic Josephson current must explain both the Kosterlitz-Thouless-type of su- 
perconductivity in the planes pHl(i|] and the effective coupling of these planes on a microscopic basis. In previous 
theoretical studies of the intrinsic Josephson effect, the supercurrent was not investigated on such a rigorous micro- 
scopic basis. Rather, it was studied using macroscopic electrodynamics employing, for example, empirically derived 
coupled sine-Gordon equations for stacked Josephson junctions [ pr| , ^2| . 

In this paper, we present numerical evidence based on Quantum-Monte-Carlo (QMC) simulations for the existence 
of the supercurrent and the intrinsic Josephson effect on a purely microscopic basis. The simulated model consists 
of a stack of coupled Hubbard planes for which, in analogy to the experimental situation, each individual plane is 
modelled by a two-dimensional (2-D) short coherence-length superconductor i.e. the 2-D attractive Hubbard 

model. 

A simplified mean-field (BCS) description of this model similar to an earlier version by Tanaka fro] is discussed in 
Sec. 0. This BCS description is similar in spirit to the conventional Josephson description in that it introduces, in 
analogy to the standard treatment, the phase dependence of the BCS order parameter (which is constant in a plane) 
perpendicular to the planes. Despite its limitations, this simple mean-field solution illustrates some of the physics: in 
particular, the dependence of the supercurrent on the microscopic parameters, i. e. the inter-plane hopping and the 
on-site interaction U. However, it cannot account for the Kosterlitz-Thouless nature of the superconductivity in the 
planes ||,|9|,|l3| and the corresponding power-law pairing correlations. 
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Sec. [II summarizes the QMC evidence for the intrinsic Josephson effect, including a finite-size study with respect 
to the number of coupled planes (up to 8). The influence of model parameters such as the temperature on the 
Josephson current is extracted. Of particular interest in view of the above-mentioned "Josephson-pairing" theory 
is that we find at best a mild elevation (for intermediate interaction U = — 4t p ) of the Josephson transition 
temperature compared to the 2-D Kosterlitz-Thouless temperature. Following an argument of Ferrcl (li), this can 
be understood within a weak-coupling (small U/t p ) BCS framework: in the normal state, single-electron tunneling 
through the barrier between the two planes lowers the energy of the system. A part of this self-energy is lost when 
the planes become superconducting because of the gap in the quasiparticle spectrum. The lost self-energy is, however, 
restored by the tunneling of Cooper pairs and the resulting Josephson coupling energy. In the Strong, Clarke and 
Anderson scenario the self-energy loss due to single-particle tunneling is assumed to be suppressed and only the energy 



lowering due to Copper pair tunneling is effective. Finally Sec. IV summarizes the results 



II. DESCRIPTION OF THE MODEL AND MEAN-FIELD EVALUATION 

The superconducting CuO^ planes in x-y-direction (Fig. [j]) are simulated by a 2-D Hubbard model with attractive 
interaction, which is subjected to periodic boundary conditions: 

H P = -t P ]T (ct aCj> + h.c.) 



<«,i>n ,a 



c i,T c L c u -A*y]ct ^ . (i) 



The summation < i,j >|| is taken over nearest neighbors in the planes and t p , U > and /i denote the transfer energy 
in the plane, the attractive interaction and chemical potential, respectively. The planes are coupled by a hopping 
term in the perpendicular (z-) direction: 

H± = -tz E c Lcj.a + h.c. (2) 

<i,j>±,a 

In the z-direction twisted boundary conditions are used. Fig. |l| displays the geometry used in the QMC simulations 
in which the Hubbard planes are stacked into a torus and threaded by a magnetic held. The magnetic held B in 
x-direction is confined to the center of the circle. Since the wave function must be single valued, a particle going once 

around the flux line acquires a phase exp ( 2 ^ ^ , where $o = ~ is the elementary flux quantum |nj and $ is the 
threaded flux. Thus, the fermionic operators are subject to the boundary conditions: 

/2tt«$\ , , 

C i+N z g z ,„ = ex P I -3— 1 %a (3) 

°t+N x S xW ,<T = C i>- ( 4 ) 

Here e x .y,z are the lattice vectors of unit length and N z and N x are the linear lengths in the z- and the ^-directions, 
respectively. Through a canonical transformation, the total Hamiltonian H = H p + H± may be written as 

H = -tp ( c *W + h - c -) 

<i,3>l\ ,a 

-*- £ ( c - c - ex K^) +kc - 

i i.a 

The fermionic operators now satisfy periodic boundary conditions in both the x- and z- directions. This model can be 
interpreted as a torus (shown in Fig. If) built up from stacks of coupled Hubbard x-y-planes and threaded by a i?-field 
in the x-direction. As shown in Sec. 01, this geometrical arrangement is directly accessible to QMC simulations. 

The influence of the B-field on the Josephson dc-current is described by the gauge-invariant form of the Josephson 
equation pljll: 
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j = sm(e 2 -9 1 -^ JiA-dl)^. (6) 

Here 82 — 8\ denotes the phase difference at zero magnetic field of the two planes considered and A gives the vector 
potential induced by B. 

In our microscopic model the bare interplane hopping t z in Eq. || is real, (i. e. the bare tunneling matrix elements). 
In this case 82 — 0± = 0, and the sinusoidal dependence of the tunneling current stems entirely from the line integral 
J A ■ dl = $, i. e. the magnetic flux: 

3 = Shl ( 4 ^) ■ (?) 

In order to verify this intrinsic Josephson behavior, the tunneling current perpendicular to the planes (in z-dircction, 
Fig. [|) can be calculated, using the relation [S3 

<Mff>- 

This current should display a sinusoidal dependence on the penetrating flux with period |. 

Let us illustrate the idea by considering first a simplified BCS treatment pS] , which is based on a Hartree evaluation 
of the free energy F and the current 

/dH\ dF 

] = \al) K WW7Y (9) 



Here, the interaction part 



is simplified to the usual BCS form 



H iat = -Uj24jC iA 4 d c id (10) 



^nt«-E( C HA + h.C.^), (11) 

where A = |Z7| (cifC^) denotes the pair potential. For an attractive on-site potential, the Cooper pairs have simple 
s-wave symmetry The usual phase dependence is introduced by writing the pair potential A m = Ae lm ^ for the 
m-th plane. 

Fourier-transforming in all spatial directions, this BCS Hamiltonian then can be written: 

H = E < k ) c L^ - E A ( fc )[4Ai + c -^ c *t] + E httt' ( 12 ) 

ka k k ' ' 

where k = (k Xl k y , k z ) is the 3-D wavevector and 

e(k) — —2t p [cos(k x a) + cos(fcj,a)] — 2t z cos [ k z a + ( ] ] (13) 

gives the kinetic energy of the electrons. This Hamiltonian, which is bilinear in its fermionic operators, can be 
diagonalized using a standard Bogoliubov transformation: 

h = E>*(«! T «fcT + 4 iQfei ) -E[^- ^} + E TP ( 14 ) 

where ik = ^[e + (k) + £-(&)], Ek = \J(\ + A 2 and Ek = Ef. + \ [t+{k) — e~(k)] is the BCS quasiparticle-energy, and 
e + (fc) = 2t p (cos(k x a) + cos(k y a)) + 2t z cos ^k z a + (j^~~^~^J^J ~ ^ 

K <I 

N~$ 



e_(/c) = 2t p (cos(k x a) + cos(k y a)) + 2t z cos [ k z a — [ I I — /i. (15) 
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The free energy can be obtained analytically from H, yielding: 



NA 2 

f = — ~ " 2fc B rE 1 °g( 1 + e ^ Ek )- ( 16 ) 

k k 

The tunneling current can be obtained by numerically differentiating F(&). The parameter A is extracted by 
minimizing the free energy, resulting in the BCS-type gap equation: 

A typical result of this simplified analysis is shown in Fig. [| which illustrates the dependence of the mean-field 
Josephson current on the microscopic parameters, t z and U (here and in what follows all energies are measured in units 
of t p ). In Fig. |2| the current j($/&o) is plotted against the magnetic flux 3>/$o f° r t z = 0.1t p at half filling and for 
various values of U. The current shows the expected sinusoidal ^-dependence with period ^. The maximum current 
decreases with increasing on-site interaction, because the effective mass of a Cooper pair and thus its localization 
tendency in a plane scales roughly with U |l5|]. If the interaction drops below a critical value, the supercurrent 
vanishes above a critical value of <&. This behavior can be understood by looking at the energy of the Bogoliubov 
quasi-particles, which can be written as 

2?r $ 



E k = ^/i 2 k + A 2 (k) - 2t z sm(k z a)sm {jj-^J ■ ( 18 ) 

If the pair potential A is smaller than the hopping energy t z , the energy Ek becomes negative for large $ and, 
therefore, a quasi-particle is excited, which corresponds to the destruction of a cooper pair. One should notice that 
for $ = 0, i. e. without an applied £?-field, Ek is positive for all finite values of A, and BCS superconductivity is 
stable for all finite values of U. Similarly, one can extract the influence of the filling (n) on the supercurrent As 
expected, scales with the number of electrons in the system. 

What goes wrong in this mean-field approach? In 2-D systems, fluctuations drive the critical temperature Tc for 



the onset of off-diagonal long-range order (ODLRO) down to zero \22\. There exists, however, a Kosterlitz-Thouless 
transition at a temperature Tkt, below which there is quasi-long-range-order (power-law pair correlations) 
Due to experimental evidence for intra- and inter-unit-cell Josephson junctions in YBaCuO single crystals |3], 
is generally believed that the high- Tc superconductivity is intrinsically 2-D in CuO^ bilayers, coupled together by 
Josephson currents along the c-axis (our z-axis) direction. There is also clear experimental evidence from the short 
coherence length (a few lattice parameters) that the 2-D superconductivity is different from the conventional BCS 
one: the typical radius ( w coherence length) of a Cooper-pair is comparable to the distance between pairs, and in 
this sense the HTSC are in an intermediate regime between the usual BCS weak-coupling regime and a regime in 
which pre-existing pairs Bose condense |^,^6| . QMC simulations |||| have correctly established the short-coherence 
length Kosterlitz-Thouless type of superconductivity for the 2-D attractive Hubbard model. It is therefore natural, 
to also apply this numerically rigorous QMC procedure to search for numerical evidence for Josephson tunneling in a 
stack of coupled CuC*2 planes. 



III. QMC RESULTS 

The Josephson current is characterized by the sinusoidal dependence j($) = Jo sin (iiT-^j. In our actual QMC 
simulation (see below) , we detect also current contributions of the form: 

/ $ 

j s ($) = J s sin 2tt— 
\ 

which are due to single-particle tunneling. This is a "finite-size" effect, with limjVj,— voo J s = 0, which has to be 
separated from the physical effect we are after. 
The current 
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perpendicular to the planes can be calculated using a finite-temperature QMC technique |2l]] for the geometry of Fig. 
|l|. The expectation value in equation (|l9| ) can directly be expressed in Green's functions Gij, which are accessible to 
standard QMC routines, i.e. 



-t z 



-it. 



<i,j>z& 

2tt 



. N z 
>, 



( n * 






-2tti — 


±) 


- h.c. 


V $0 


nJ 





a<F 



El $ 1 
cLc^vr exp(— 2ni ) — h.c. 

<i j> 2 ,ct 



(20) 



where $ = ^4 2 L with L — N z a and a is the lattice parameter. Using 

( C i<j C j<?) ~ ~ ( C j<7 C ia) ~ ^ij ~ G 'jit 

the current becomes: 



2tt \ - 

0) = -^ E 



a$ 



<i,3>z<J 



$ 1 \ / $ 1 

G% exp ( 27rf — — ) - GJ, exp ( -2**— — 



(21) 



In the first simulation we used a 4 x 4 x 4 lattice with periodic boundary conditions in all three spatial directions. 
Because of the limited system size, the results display considerable finite-size effects. In order to separate these finite- 
size effects from the physical effect we are after, we have calculated j($>) in the non-interacting (U = in eq. (0)) 
tight-binding model for two different sizes (Fig. ||). The results show, that the finite-size effects are also sinusoidal, 
but with period 1, and, indeed (see the argumentation below), they are related to single- particle tunneling. Thus, 
if a Josephson current appears in the system, j(<&) should behave like a superposition of two sinusoidal curves with 
period \ and 1. It should, therefore, be of the following additive form 



= J 3 sin ( 4tt— 



J, sin ( 2tt — 
$ 



(22) 



with the first term referring to the Josephson supercurrent and the second to the single-particle finite-size effects. The 
parameters Jj and J s are the corresponding amplitudes. 

The Hubbard model parameters used in the following discussion are: U = — 4t p , (n) — 0.75, t z — 0.1£ p , for which 
the 2-D attractive Hubbard planes are known to show Kosterlitz-Thouless type of superconductivity for /3 = 8 - 10 
@|- 

As seen from Fig. |4|, which contains the QMC results (diamonds plus errorbars), the curve ( |22j ) fits the calculated 
points very well. Thus, we have a first evidence that our model indeed describes a Josephson junction. 

By increasing the number of planes from 4 to 8, the amplitude, J s , of the single-particle finite-size effects can be 
reduced (Fig. ||), and the sinusoidal curve with period ^ can be seen more clearly. The interpretation of j s ( < E > ) = 

J s sin (^2ir-^j as a finite-size contribution, that is related to the single-particle tunneling, is based on its scaling 

behavior as a function of the number of N z of planes. The reason for this finite-size behavior of j s can simply be 
understood: Enhancing the number of coupled planes, N z , the phase of the single-particle hoppings in eq. (||) is more 
effectively averaged out to zero than for smaller N z , and J s should scale to zero. 

On the other hand, Jj, the Josephson amplitude, should scale with the number of planes N z , a relation, which is 
clearly obeyed by our QMC data in Figs. [| and |^. In Fig. ^ and [|, ^-j ( < I ) ) is plotted against <|- for 4x4x4 and 
4x4x8 systems, respectively. The lines indicate the interpolated curve. As expected, the magnitudes of jj-Jj in 
Figs. ^ and || are equal {Jj — 0.0079 is used in both figures). 

The 2-D attractive Hubbard model for U = —4 and (n) = 0.75 shows a Kosterlitz-Thouless transition at j3 c t p « 10 
from the superconducting to a normal phase This transition should be reflected in coupled Hubbard planes, and 

the Josephson current should break down at some (3 — [3 C . The temperature dependence of j allows an estimation of 
(3 C . In Figs. and||, is plotted for (3 = 2, 6 and 7 respectively. Fig. ^ gives the corresponding plot for (3 = 8. For 
(3 = 2, is sinusoidal with period 1, i. e. no Josephson current appears, and the system is not superconducting. For 
(3 = 6 and (3 = 7, Jj becomes finite, but is still small compared with the value of Jj for (3 = 8. The J values for (3 = 10 
(not plotted) have been found to be essentially unchanged compared to (3 = 8. This growing current (with saturation 
around (3 ss 8) is due to interlayer tunneling of local pairing of electrons. For temperatures (3 < 8, exponential not 
power-law pairing correlations are present in the 2-D Hubbard model. Accordingly, preliminary results indicate that 
for (3 < 8 the current does not scale with the number of sites in a plane. A rough estimate for the Josephson transition 
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temperature is, therefore, = 8, at best a mild elevation compared to the single plane. Within BCS Ferrel showed 
that the shift in ground-state energy due to tunneling is the same in the normal and superconducting stats fl4|| . The 
physical significance of this result is that the ground-state energy of the junctions in its superconducting state is 
neither favored nor disfavored, relative to the normal-state energy by the presence of the tunneling barrier. Therefore, 
in this BCS limit, we do not expect any elevation in the transition temperature. 

IV. CONCLUSIONS 

In this paper, we have presented a model for intrinsic Josephson couplings on a microscopic length-scale. It is 
based on coupled Hubbard planes with attractive interactions, which are threaded by a local magnetic field. The field 
dependence of the current perpendicular to the planes is calculated using numerically rigorous Quantum-Monte- 
Carlo simulations. To extract some of the physical ideas, in particular the dependence on the microscopic parameters 
U/tp and t z /t p , we first reviewed the mean-field BCS approximation. The maximum current j max increases with 
decreasing on-site interaction U and increasing hopping energy t z . At small values of U , for which the pair potential 
A is smaller than t x , the superconductivity is unstable for finite values of the magnetic flux $. In mean-field, the 
Hubbard model reduces to an effective BCS type Hamiltonian, which is not capable of correctly describing salient 
features of the 2-D superconductivity (short coherence length, Kosterlitz-Thouless transition) in the CuOi planes of 
the HTSC's. 

To overcome these shortcomings and to present a rigorous proof of the intrinsic Josephson couplings, we have 
extracted from the, in principle (apart from controlled statistical error |H]), exact Quantum-Monte-Carlo sim- 
ulations. The results show a sinusoidal behavior of A careful discussion of the finite-size scaling of the current 
clearly reveals that it is the behavior of a Josephson supercurrent. The Josephson current drops to zero for j3 < 6, at 
U = —4 and a filling (n) — 0.75. Due to finite-size effects, this transition is smeared out and, therefore, we can only 
approximatively determine the transition temperature from superconducting to the normal state at /? = 8. These 
facts give clear evidence that the attractive Hubbard model is capable of describing the Josephson effect in HTSC's 
on microscopic length scale. Careful studies of the parameter dependence of the maximum supercurrent are under 
way. These studies should enable us to systematically extract the influence of material properties on the Josephson 
current, and thus should be useful to optimize the intrinsic Josephson effect for specific applications. 
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FIG. 1. Schematic view of the used geometry. The coupled superconducting planes are arranged in a torus threaded with a 
the local magnetic field B. 
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FIG. 4. Current perpendicular to the plane for the 4x4 (planar size) x4 (number of planes) system QMC-results: diamonds 
plus errorbars, the dotted line shows the interpolated curve j($). 
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FIG. 5. Current perpendicular to the plane for the 4 x 4 x 8 (8 coupled planes in Fig. [IJ) system. The dashed line shows the 
U = behavior. 
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FIG. 6. Current perpendicular to the plane for (3t p = 2. 
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FIG. 7. Current perpendicular to the plane for the same parameters as in Fig. |6|, however for a three times lower temperature, 
i. e. (3t p = 6. 
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(3 = 7 4x4x4 U = -4 (n) « 0.75 
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FIG. 8. Current perpendicular to the plane for (3t p = 7. 
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